Phase transition induced hydrodynamic instability and Langmuir-Blodgett Deposition 
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We propose a model to understand periodic oscillations relevant to the origin of mesoscopic channels 
formed during a Langmuir-Blodgett deposition observed in recent experiments {M. Gleiche, L.F. 
Chi, and H. Fuchs, Nature 403, 173 (2000)}. We numerically study one-dimensional flow of a van 
der Waals fluid near its discontinuous liquid-gas transition and find that steady-state flow becomes 
unstable in the vicinity of the phase transition. Instabilities leading to complex periodic density- 
oscillations are demonstrated at some suitably chosen sets of parameters. 
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Studies of pattern formation in non-equilibrium sys- 
tems have attracted considerable attention due to their 
relevance in understanding fundamental physics and the 
potential for applications in a wide range of emerging 
technologies. The spontaneously formed patterns can 
be converted into functional structures when the pattern 
forming mechanism is properly controlled and designed. 
Typically, this class of pattern formation involves sys- 
tems flowing at conditions around which the properties 
undergo dramatic changes. The patterns are a result of 
the occurrence and propagation of a phase transition or 
chemical reaction, the Belousov-Zhabotinsky reaction |l]] 
and capillary flow of molten polymer || being two ex- 
amples. One of the key issues to be understood in these 
systems is the nontrivial coupling between the system 
specific 'phase transition' and the hydrodynamics. 

Our work has been motivated by recent experiments 
that have demonstrated the formation of macroscopic 
arrays of submicron-sized channels during a Langmuir- 
Blodgett deposition S. Mechanisms based on stick- 
slip conditions have been proposed to describe the phe- 
nomenon H|J. However, in this Letter, we propose a dif- 
ferent scenario in which density oscillations occur in the 
flow of a single-component fluid at conditions tuned ap- 
propriately near its discontinuous liquid-gas phase tran- 
sition. The system fails to support steady-state flow at 
these conditions due to the absence of a mechanically 
stable uniform state in a certain density interval. The 
instability leads to periodic oscillations in density or al- 
ternating appearance of liquid and gaseous phases in the 
system. The oscillations can be potentially useful in pro- 
ducing large arrays of periodic structures and we propose 
this as a possible mechanism for the experiments || men- 
tioned above. The problem of fluid flow in the vicinity of 
a phase transition has been widely studied M . A closely 
related subject on propagation of the phase transition 
front has also been investigated in detail Q. However, 
we are not aware of any previous work on unstable one- 
dimensional flow of a single-component fluid that leads 



to periodic density-oscillations of this nature. A realiza- 
tion of such a system is a fluid flowing in a narrow chan- 
nel, or higher dimensional flow in which spatial variations 
of the system properties are suppressed in the direction 
transverse to the fluid velocity. The flow of surfactant 
molecules in a typical Langmuir-Blodgett deposition is 
an example in two-dimensions 

Our analysis begins with the standard hydrodynam- 
ics at constant temperature consisting of the continu- 
ity equation and the Navier-Stokes equation for a one- 
dimensional flow of compressible fluid as follows: 
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where p is the density, v is the velocity and r\ is the 
viscosity. The main feature of the system is contained in 
the quantity P, the local pressure of the inhomogeneous 
fluid. It is expressed as M 
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which is also the equation of state of the fluid. Here, the 
quantity ^[p] denotes the Helmholtz free energy density 
for an inhomogeneous fluid given by, 




ip(p) } dx, 
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where ip(p) is the Helmholtz free energy density of the ho- 
mogeneous system and k is a phenomenological constant 
associated to the gradient term with a meaning that is 
specific to the particular system chosen. For a van der 
Waals (vdW) fluid, the Helmholtz free energy density of 
a uniform system is given by 



ip(p) = SpT\n 
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inlet at x = 

v(0,t) = Vi 
p(0, t) = pi 



outlet at x 
v(L,t) = v Q 
p(L,t) 



L 



FIG. 1. Schematics of the flow configurations. The quanti- 
ties Vi and pi are, respectively, the inlet velocity and density 
while the outlet velocity is denoted by v . 



where T is the temperature and K is proportional to the 
surface tension of the liquid-gas interface. The reason 
for choosing a vdW fluid for our investigation is that it is 
one of the simplest descriptions that includes a discontin- 
uous liquid-gas transition and a region in the P — p plane 
where the uniform fluid is unstable (negative compress- 
ibility) below the critical temperature T c = 1. We note 
that the quantities P, p and T in Eqs. (||) and (^) have 
been normalized by their respective values at the critical 
point for convenience. We also remark that the fluid is 
maintained at thermal equilibrium at all times, that is, 
Eq. (||) always holds. It is the mechanical stability that 
is being investigated here. 

A numerical algorithm has been developed to solve 
Eqs. (0) and (§) using the forward-time-centered-space 
approach described in Ref. Jig] . We will restrict our- 
selves to a special case depicted in Fig. [TJ: a flow in the 
region x E [0, L] subject to the following set of boundary 
conditions 



v(x < 0, t) = Vq, 
v(x > L, t) = v L , 
p(x < 0, t) = po, 
p(x >L,t) = p{L, t). 
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The fluid velocity has been arbitrarily chosen so that it 
flows from left to right. The positions x — and x = L 
correspond to, respectively, the inlet and the outlet. The 
quantities Vq and vl are, respectively, the velocities at 
the inlet and the outlet. The density at the inlet is kept 
at po- This set of boundary conditions describes a situ- 
ation where both the fluid density and velocity are kept 
fixed at the inlet, and the fluid density at the outlet is 
allowed to vary with time while keeping the velocity con- 
stant. Under this special circumstance, the instability of 
the flow is necessarily reflected in the time dependence 
of the density at the outlet and we systematically ana- 
lyze p(L, t). We remark that there are a total of six 
parameters involved, namely, vq, Vl, po, T, rj and k. In 
general, the solutions to Eqs. ([!]) and (|2|) are complex 
and possess many degrees of freedom. Here, we will fo- 
cus our discussions on the results for a restricted sets of 
parameters. 

When the density at the inlet po and the velocity at 
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FIG. 2. The comparison between the computed instability 
boundary (blue) and that estimated from the spinodal bound- 
ary (red) for p = 0.8, v L = 0.1, r) = 0.01 and k = 0.0001. 
Also shown in green is the boundary between the simple and 
complex oscillatory behavior in the unstable region. 

the outlet vl are kept fixed, the behavior of p(L, t) as 
i>o varies can be understood intuitively as follows: if one 
insists on a steady-state flow, the density at the outlet 
p(L, t) has to be 



p(L, t) 



vqPo 
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Since Eq. (||), the equation of state of the fluid, allows 
for a region where the compressibility is negative (spin- 
odal region), mechanical instability is thus inevitable if 
vqPo/vl falls within the spinodal region. A rough es- 
timate of the stability boundary, assuming the system 
is uniform, is given by vl's for which p(L, t) coincides 
with the spinodal boundary (i.e., state of the fluid where 
compressibility diverges). Our numerical results indeed 
show that the steady state solution fails to be stable in 
some regions of the parameter space, restricted to the 
v o — T plane in this work. The instability region in the 
vq — T plane has been identified numerically. Figure § 
compares the computed stability boundary (blue) and 
that estimated from the spinodal boundary (red). The 
discrepancy between the estimates and the numerical re- 
sults is a consequence of the shift in the spinodal region 
due to the flow-induced density gradient and is evident 
from Eq. (@). 

Our results further show that the instability leads to 
oscillatory behavior in p(L, t), the density at the out- 
let. The density-oscillations have been classified into two 
classes: a 'simple' oscillation in which p(L, t) oscillates 
between two densities, and a more 'complex' oscillation 
where p(L, t) cycles among three or more densities. The 
left panels of Fig. || illustrate the behavior of a sim- 
ple density-oscillation whereas the right panels show a 
case of complex oscillations. Figure ||(a) displays plots 
of p(L, t), showing the time evolution of the density at 
the outlet. The plots of spectral density \p(L, ui)\ 2 are 
shown in Fig. [|(b). A primary peak accompanied by 
its weaker harmonics is shown on the left panel; a more 
complicated structure exists on the right panel. Further- 
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FIG. 3. Plots of (a) p{L, t), (b) \p(L, u)\ 2 and (c) P(L, p) 
for the case of simple oscillation at vo =0.142 on the left panel, 
and for the case of complex oscillations at no = 0.136 on the 
right. The other parameters are set to vl = 0.1, po = 0.8, 
T = 0.987, t] = 0.01, and « = 0.0001. 



more, the number of densities among which p(L, t) cy- 
cles can be picked out by evaluating the density distri- 
bution function P(L, p) = TV -1 j S[p(L, t) — p]dt where 
Af = J P(L, p)dp is the normalization constant. Figure 
H(c) displays P(L, p) with two peaks on the left and three 
on the right. The green curve in Fig. |^ is the boundary in 
the Vq — T plane separating the region of simple and com- 
plex oscillations. Simple oscillations occur in the region 
between the blue and the green curves, while complex 
oscillations are found under the green curve. 

Figure f|(a) is a density plot of \p(L, oj)\ 2 versus vq and 
u>, and Fig. |](b) is a density plot of P(L, p) as a function 
of vq and p, showing the dependence of these functions 
on the inlet velocity vq, while keeping other parameters 
fixed. Higher intensity is represented by darker colors 
in the density plots, highlighting the peaks of the func- 
tions. The appearance of peaks at nonzero frequency 
in \p(L, uj)\ 2 , and the multi-peak structures in P{L, p) 
reflect the unstable region in v [v G (0.119,0.146)] 
at T — 0.987. The region of complex oscillation at 
T = 0.987 with v G (0.123,0.139) can be identified as 
the deviation from the simple structure of a single prin- 
cipal peak in |/5(L, u>)\ 2 [the trend of the principal peak 
is shown by the pink curve in Fig. ||(a)], and the onset of 
additional peaks in P(L, p) other than the two expected 
for simple oscillations. 

The scenario described by Eqs. (Q) and (g) with bound- 
ary conditions given by Eqs. (||)-(g) can be related to the 
flow of a Langmuir monolayer (a single molecular layer 
of insoluble surfactant molecules spread at the air/water 
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FIG. 4. Density plots of (a) \p(L, lo)\ 2 as a function of lo 
and vo, and (b) P(L, p) as a function of p and wo. The pink 
curve in (a) is a guide-to-eye for the trend of the principal peak 
in the region of simple oscillations. Deviation from this curve 
indicates the region of complex oscillation vo G (0.123, 0.139). 
This region can be clearly identified in (b) by the onset of 
additional peaks. The parameters used are po = 0.8, vl = 0.1, 
T = 0.987, T) = 0.01 and n = 0.0001. 



interface) during a Langmuir-Blodgett deposition. In a 
Langmuir-Blodgett deposition, a substrate is withdrawn 
from or dipped into the water transferring the surfactant 
molecules onto the substrate at the contact line. The 
surfactant molecules are removed at the substrate with- 
drawal speed at the contact line. The boundary condi- 
tions at the outlet, Eqs. (Q) and (||), are explicitly real- 
ized at the contact line if the transfer ratio is unity at 
all times. Although the correspondence between the in- 
let boundary conditions and the actual experiments is 
not completely clear, Eqs. @ and (||) are reasonable 
for an appropriately chosen L. A sensible estimate of 
the inlet-outlet separation L for the Langmuir-Blodgett 
deposition would be the meniscus height at the contact 
line, which is ~ 10~ 3 m. Hence, the specific choice of 
boundary conditions, Eqs. (§|)-(§|), can be realistic. De- 
spite the attempt to associate the proposed scenario to 
the actual experiment, the model is far from a complete 
description of the system. Many details relevant to the 
Langmuir monolayer such as long-range dipolar interac- 
tions between the surfactant molecules, surfactant-water 
interactions ||, water flow near the contact line j7|,pT|, 
elasticity of the monolayer Jl2] |, are not taken into ac- 
count. Nevertheless, the key emphasis is to point out 
that a discontinuous phase transition has been shown to 
lead to simple period density-oscillations which provide 
a possible mechanism for the formation of macroscopic 
arrays of mesoscopic structures. 

In summary, we have demonstrated numerically that a 
one-dimensional flow of a van der Waals fluid near its dis- 
continuous liquid-gas transition exhibits a hydrodynamic 
instability leading to periodic, classified as simple and 
complex, oscillatory behavior. We highlight the key point 
of our analysis: the instability in this simple system arises 
because of the inability to support a steady-state solution 
when enforcing the continuity equation in the vicinity of 
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a discontinuous transition. Hence, this oscillatory be- 
havior is not restricted to van der Waals fluids and could 
be observed in other kinds of systems (with well pre- 
scribed equation of state) that possess similar features. 
We have attempted to relate this one-dimensional model 
to Langmuir-Blodgett deposition where the boundary 
conditions can be physically realized. This simple model 
provides us with a route to periodic oscillatory behavior 
that has been observed S. We note that the region in 
the parameter space where simple oscillations occur is 
limited. This is a possible explanation for the difficulty 
in finding the periodic oscillations and their sensitivity to 
experimental conditions, such as pressure and tempera- 
ture. Estimates for the period and the interval of each of 
the phases, within the framework of this model, require 
a knowledge of the continuation of the equation of state 
into the metastable regime and the spinodals. 
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